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STABILITY OF OAS TIRING- WHEELS FOR AIRCRAFT LANDING- GEARS 



By Arthur Kantrowitz 



SUMMARY 



A theoretical study has "been made of the shimmy of 
castering wheels. The theory is based on the discovery of 
a phenomenon called "kinematic shimmy." Experimental 
checks, using a model having low-pressure tires, are re- 
ported and the applicability of the results to full scale 
is discussed. Theoretical methods of estimating the spin- 
dle viscous damping and spindle solid friction necessary to 
avoid shimmy are given. A method of avoiding shimmy - lat- 
eral freedon - is introduced. 



INTRODUCTION 



In many installations of castering rubber-tired wheels 
there is a tendency for the wheel to oscillate violently 
about the spindle axis* This phenomenon, popularly called 
"shimmy," has occurred in some airplane tail wheels and 
has-been corrected in two ways: First, by the application 
of friction in the spindles of the tail_wheels; and, sec- 
ond, by locking the wheels while taxying at high speeds. 
Shimmy is common with the large wheels used as nose wheels 
in tricycle landing gears and, since it is impossible to 
lock the wheels, friction in the nose-wheel spindle has 
been the sole means of correction. Because the nose wheel 
is larger than the conventional tail wheel and usually car- 
ries a greater load, the larger amounts of spindle friction 
necessary to prevent shimmy are objectionable. 

Several experimental investigations undertaken to find 
methods of avoiding shimmy showed that shimmy could be pre- 
vented in small tail wheels by friction in the spindle. 
The most thorough work on the general subject of stability 
of castering wheels is that of Becker, Fromm, and Maruhn 
(reference l) on shimmy in automobiles. Because they do 
not include the effect of lateral tire deflection, it is be- 
lieved that their results would, in general, not be appli- 
cable to airplanes. 

The present paper presents a theoretical and experi- 
mental study of the problem of the stability of castering 
wheels for airplane landing gears. On the basis of simpli- 
fied assumptions induced from experimental observations, a 
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theory of the phenomenon is presented. The theory is then 
compared quantitatively with the results of model experi- 
ments. 

THEORY OF THE STABILITY OF CASTER ING WHEELS 
Kinematic (or Static) Shimmy 



Some preliminary experimental results on shimmy were 
obtained "by the N.A.C.A. with the aid of the "belt-machine 
apparatus shown in figure 1, This machine consists of a 
continuous fabric belt mounted on two rotating drums and 
driven "by a variable-speed electric motor. Provision is 
made for rolling a castering wheel up to about 6 inches 
in diameter on the "belt in such a way that it is free to 
move vertically "but not horizontally. 

On this "belt machine the following phenomenon (see 
fig. l) was discovered while pushing the holt very slowly 
"by hand. With the wheel set at an angle with the "belt as 
in (a) and the "belt pushed slowly, the "bottom of the tire 
would deflect laterally as is shown in ("b). When the "belt 
was pushed farther, the wheel straightened out gradually 
as is shown in (c). However, the "bottom of the tire tfould 
then still "be deflected and the wheel would continue to 
turn as in (d). The wheel would. thus finally overshoot, 
as shown in (e) and (f). The process would then he re- 
peated in the opposite direction. 

Figure 2 is a photostatic record of the track left "by 
the "bottom of the tire on a piece of smoked metal. Two 
things will he noticed: First, that the "bottom of the 
tire did not skid; and second, that the places where the 
wheel angle is zero (indicated "by zeros on the track) cor- 
respond roughly to the places where the lateral deflection 
of the tire is a maximum. Thus the wheel angle lags the 
tiro deflection by l/4 cycle. 

It" was noticed that the oscillation could he inter- 
rupted at any point in the cycle "by interrupting the mo- 
tion of the belt without appreciably altering the phenom- 
enon. From this observation it was deduced that dynamic 
forces play no appreciable part in this oscillation. This 
type of oscillation might be called "static" shimmy but 
the more rigorous term "kinematic shimmy" is preferred. 

It was also found that the distance along the belt 
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required for one cycle did not vary much with caster angle 
or caster length. (Sec fig. 3.) Caster length was there- 
fore considered not to "be of fundamental importance in 
this type of oscillation. 

It should "be pointed out that , in order to observe 
the kinematic shimmy, lateral restraint of the spindle is 
necessary to prevent the spindle from moving laterally 
when the bottom of the tire is deflected and thus neutral- 
izing the tire deflection. This restraint is supplied by 
the dynamic reaction of the airplane when the airplane is 
moving forward rapidly but is not ordinarily present when 
the airplane is moving forward slowly. It has been ob- 
served, however, on airplanes towed slowly with two tow 
ropes so arranged as to provide some lateral restraint. 

It is observed in figure 1(a) that, when the center 
line of the wheel is at an angle 6 (see fig. 3) with the 
direction of notion, the bottom of the tire deflects. 
This situation is represented schematically in figure 4. 
It is seen that a typical point on the peripheral center 
line must have a component of motion perpendicular to the 
wheel center line if the tire is not to skid. Thus 

&A = r sin 6 ds 

(The minus sign follows from the conventions used as shown 
in fig. 3.) Since only small oscillations are to be con- 
sidered, the approximation 

Q = ~ e (i) 

ds 

may' be substituted. 

The effect of tire deflection on 6 will now be con- 
sidered. It will be assumed, for the purposes of rough 
calculation, that as illustrated in figure 5 the projec- 
tion of the peripheral center line on the ground is a cir- 
cular arc intersecting the wheel central plane at 'the ex- 
tremities of the projection of the tire diameter. (See 
fig. 1(d).) Thus, in figure 5, r is the tire radius. 
(It will be assumed for the time being that the castor 
length and caster angle arc zero.) Now if the tire is de- 
flected in the forn of a circular arc, then the condition 
that the torque about the spindle axis be zero is that the 
strain bo symmetrical about the projection of the wheel 
axle on the ground. Clearly, if the wheel is displaced, 
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it will be turned about the spindle axis by the asymmet- 
rical elastic forces until, if it is allowed time to reach 
equilibrium, the symmetrical strain condition is reached. 
Thus, if the tire is deflected an amount A as in figure 
5(a); and, if the wheel rolls forward a distance ds to 
the condition shown in figure 5(b), in order- for the strain 
to remain symmetrical the wheel must turn about the spin- 
dle axis an amount d9 . From figure 5, E d8 = & s . The 

readily obtained from geometry in 



R 
r 



may 
and 



value of 
terms of 

which, if X s 
substituting for 



« r 



be 
A. 

a 



Thus R s = r 2 + (R - X) 3 ; from 
it is seen that R = r 2 /2X. Then 



R. 



dG_ 

ds 



(2) 



If the caster length is finite, the strain will not 
be symmetrical about the axle as was assumed here but will 
be symmetrical with respect to some line parallel to the 
axle but a certain fixed distance ahead. Hence, the es- 
sential elements of the geometry are unchanged and all the 
reasoning that led to equation (2) is still valid for 
this case. 

Since the phenomena represented by equations (l) and 
(2) occur simultaneously, they must be combined to get 
the total effect. Thus 



d 2 e 
ds 2 



_2_ 

„2 



0 



(3) 



This differential equation corresponds to a simple 
harmonic oscillation occurring every time the wheel moves 

distance 2tj/,J2/t 3 = tt r. 



a 



Measurements of the interval of kinematic shinny S 
have been made for three similar tires of the type illus- 
in figure 1. These tires all had radii of approx- 
so that the theoretical interval was 

Their experimental intervals were 0,65 
and 0.79 foot. This agreement is closer 
been anticipated in view of the roughness 



trat ed 
imately 2 inches 
about 0.74 foot, 
f oo t , 0.74 f oo t , 

might 



than might have 

of the assumption. It will be seen from equation (l) that 
A is l/4 cycle out of phase with G , which is in agreement 
with the information obtained from figure 2. 



In order 
s umx> t i 0 n made 



to take account of tires for which the as- 
concerning the projection of the peripheral 
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center line is not quant itatively valid, an empirical con- 
stant K 1 will be used in place of 2/r 2 in equations 
(2) and (3), thus obtaining 

M = Kl A (4) 
ds 

and 

= - K x e (5) 

ds 2 

The constant K x can be measured by observing the inter- 
val of kinematic shimmy. Where experimental values of K x 
are available, they will be used rather than the rough 
theoretical value 2/r 2 . 

Dynamic Shimmy 

In the foregoing derivation for the oscillation called 
"kinematic" shimmy it was assumed that the strain of the 
tire was always symmetrical, that is, the wheel was moving 
so slowly that any torque arising from dynamic effects in- 
volved in the oscillation would be negligible. For this 
case, from equation (4) 

\ Ki ds/ 

If, now, the wheel is assumed to be moving at a ve- 
locity such that the effect of the moment of inertia about 
the spindle axis is significant, then the strain can no 
longer be symmetrical and for small asymmetries the torque 
exerted by the tire on the spindle will be proportional to 
the amount of the asymmetry. Thus, the value in paren- 
theses will no longer be zero but it can be assumed that 
it will be proportional to the dynamic torque, 




where Ci is an appropriate constant of proportionality 
and includes the moment of inertia, about the spindle axis 
I w . If the forward velocity 7 of the wheel is constant 

d 8 6 ? a dfjB 
dt 2 ds 2 

and 
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v= A -_ o, (x - i 3£) (6) 
ds V K, ds> 



The constant C x may "be neasured by deflecting the 
"bottom of the tire a known amount, moving the wheel for- 
ward, and balancing the torque T exerted by the tire on 
the spindle so that 6 stays constant. The method of de- 
flecting the bottom of the tire a known amount will be 
described later. In this case (6) becomes 

- T - = C x X 
x w 

from which 0i may be found. It has been found that Q x 
increases with increasing caster angle. Thus, for a tire 
like the one in figure 1, C x was 71,000 radians per 
second 2 per foot for a caster length L of 0.17 inch 
(caster angle, 5°; no fork offset) and was 104,000 radians 
per second 8 per foot for caster length 0.68 inch (caster 
angle, 20°). 

In the study of kinematic shimmy it was also seen that 
the only change in X was due to the fact that a compo- 
nent of the forward motion was perpendicular to the cen- 
tral plane of the wheel. This circumstance is expressed 
by equation (l). It is, however, found that, if the spin- 
dle is clamped at 6=0 and the bottom of the tire is 
deflected, the deflection will gradually neutralize itself; 
that is, the bottom of the tire will roll under the wheel. 
Thus the asymmetrical (d6/ds = 0) strain that exists in 
this case contributes to dX/ds. The case of 6 and 
dG/ds aero is illustrated in figure 6. If the effect is 
again supposed to be proportional to the cause, there is 
obtained 

= - G 3 X (?) 

ds 



In this equation the constant 0 2 is a geometrical 
constant of the tire that can be obtained from static 
measurements. The order of magnitude of C 3 can be ob- 
tained by assuming that the periphery, of the tire inter- 
sects the extremity of the extended central plane of the 
wheel. In that case C 2 = l/r. 

If G is not zero, there will be a component of the 
forward notion contributing to dx/ds. As in equation 
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(1), this component will be -G. Adding this component 
to the part of dX/ds due to asymmetry, then (d6/ds 
still assumed zero) 

4- = - e - c 2 x (s) 

ds 

d0 

This equation expresses that for - = 0 the contri- 
bution of the asymmetrical strain to dX/ds was - C 3 X • 

Also for the symmetrical strain, in which case X - = 

0 (kinematic shimmy), there was, of course, no contribu- 
tion due to asymmetry. Assume now a linear interpolation 
"between these two limiting cases. Thus, finally, 



ds 3 \ Ki dsy 



A method of measuring G 3 is provided by equation 
(7) which, when integrated, gives 

log e v- = 0 a (s - s 0 ) 
o 

If, with the spindle clamped at 6=0, the tire is 
deflected a known amount X Q and rolled ahead a known 
distance and the new X measured, C 2 may be computed. 
It was seen earlier that C2 was of the order of magnitude 
of l/r, that is,- it would be of the order of 6 for a 2- 
inch radius tire. The constant C 2 was measured on two 
model tires under different loads and found to be 6.2 and 
3.4. Considerable variations of this constant with tire 
pressure and load have been found. 

A method of obtaining the constant known X necessa- 
ry for the measurement of Ci is provided through equa- 
tion (8). Here it is seen that, if the wheel is pushed 
along at a constant angle 6, X will increase (negative- 
ly) until 8 = - G 2 X , in which case J^A = 0 and equi- 

ds 

librium is reached. 

When the wheel is moving ahead at a finite constant 
Velocity, the phenomena represented by equations (6) and 
(9) occur simultaneously. Therefore, to get the total ef- 
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foot, combine the two equations thus obtaining 

il £l + (± + !* ^ £| + 6 = 0 (10) 

The solution for the natural modes of motion repre- 
sented "by equation (10) is 

o&i s a s s a,s . > 

B=Ae + B e + C e 3 (11) 

where the a's are the three solutions of the so-called 
"auxiliary" equation 



(_L + °£ v a ") a 5 + 1 = 0 (12) 
\K 1 0 1 /' 



One of these a's is real and negative arid corre- 
sponds to a nono scillatory convergence. The other roots 
are conjugate complex numbers and correspond to the shimmy 
under consideration. The roots will "be of the form 
a ± O) i. If the divergence a, is positive, the oscilla- 
tion will steadily increase in amplitude (while its ampli- 
tude is not large enough for skidding to occur); and, if 
a is negative, the oscillation will steadily decrease in 
amplitude and eventually disappear. The meaning of the 
quantity "divergence" may he illustrated by saying that it 
is approximately equal to the natural logarithm of the ratio 
of successive maximum amplitudes to the distance between 
them. The quantity u> is equal to 2tt times the number 
of oscillations per foot. The frequency therefore is 
u>V/2tt, The phase angle is obtained by substituting for 6 
in equation (6) the value obtained from the foregoing pro- 
cedure and solving for A. 

The divergence, 'the frequency, and the phase relations 
thus derived for typical model tire constants are plotted 
in figure 7. For small velocities (0 to 6 feet per second) 
the frequency corresponds to kinematic shimmy; it is pro- 
portional to velocity. The divergence increases rapidly, 
however, because the spindle angle lags on account of the 
moment of inertia about the spindle axis, thus allowing 
more lateral deflection than would occur in a kinematic 
shinny. On the next half cycle a larger spindle angle is 
reached and the process repeats. As the velocity is fur- 
ther increased, the lag, and hence the asymmetry of the 
strain, further increase until the strain becomes almost 
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entirely asymmetrical. For this condition, ( — ) — « A. 

^ds' Ki 

Then the restoring torque on the spindle is approximately 
proportional to A. (See equation (6).) However, X 
(measured negatively) will still lag somewhat "behind .8 
"because, after the wheel is turned through a given angle, 
a certain forward distance is required for equilibrium 
tire deflection to "be reached. Thus, the restoring force 
will again lag the displacement. As the velocity increases 
in the high-velocity range, the frequency stays nearly con- 
stant (see fig. 7) and the distance corresponding to a sin- 
gle oscillation increases. Hence, this constant lag be- 
comes a smaller part of the cycle and the divergence de- 
creases at high velocities. 

It will be appreciated that the foregoing theory con- 
siders only .the fundamental phenomena taking place in 
shimmy. There are, as would be expected, other phenomena 
occurring s iraultaneously which have been neglected. Some 
of the. more important of the phenomena neglected are: 

1. Miscellaneous strains (other than lateral tire 
deflection) occurring in the tire. A' rubber tire being an 
elastic body will distort in many complicated ways while 
shimmying. In particular, there will be a twist in the 
tire due to the transmission of torque from the gro^und to 
the wheel. 

2. Two effects will cause the stiffness constants of 
the tire to change with speed. First, centrifugal force 
on the rubber will make the tire effectively stiffer at 
high speeds. Second, much of the energy used to deflect 
the tire will go into compressing the. air. The compressi- 
bility of the air will change with the speed of compres- 
sion owing to the different amounts of heat being trans- 
ferred. from it. 

3. There will be a gyrostatic torque aoout the spin- 
dle axis caused by the interaction of the rotation of the 
wheel on the axle and the effective rotation of part of 
the tire about a longitudinal axis on account of the later- 
al tire deflection. This torque will later be shown to 
have a noticeable effect on the results. 

The inclusion of items 1 and 2 in the theory would 
obviously be very difficult. It is therefore necessary to 
resort to experiment to determine whether the present the- 
ory gives an adequate description of the phenomena. If 
so, the omission of these and any other items will be jus- 
tified. 
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An experimental check on the theory was obtained "by 
measuring the divergence and frequency of the shimmy on 
the "belt machine at two caster angles and at a series of 
velocities. These measurements were made "by placing a 
lighted flashlight bulb on a 6-inch sting ahead of a model 
castering wheel with a ball-bearing spindle and then tak- 
ing high-speed moving pictures of the flashlight bulb with 
the wheel free. The photographs were made with time re- 
cordings on the film, and the belt carried an object that 
interrupted light from a fixed flashlight bulb and thus 
recorded the belt speed on the film. 

The divergence and frequency of the shimmy were ob- 
tained by measuring the displacements and times corre- 
sponding to successive maximum amplitudes (while the ampli- 
tude was still small enough to make all the assumptions 
valid). The results are plotted in figure 8. 

In order to compare these results with the theory, 
the constants C x , C2 , and Ki were measured on the same 

tire at the two caster angles by the previously described 

methods. It was found that G 2 = 6.2 feet" 1 ; K x = 62.5 

feet"* 2 ; that for 5° caster angle, C x = 71,100 feet" 1 - 

second^ 3 , and for 20° caster angle, Oi = 104,000 feet" 1 - 

second"" 2 . The roots of equation (12) were then found and 
the divergence and frequency of the shimmy computed for a 
series of velocities and at caster angles of 5° and 20°. 
These results are plotted in figure 9 and, for purposes of 
comparison, the experimental curves are also reproduced. 

The agreement between theory and experiment is consid- 
ered satisfactory as regards qualitative results. It will 
be noticed, however, that the theoretical values of the di- 
vergence are decidedly too large at high velocities, say 
25 feet per second. It should be noted that 25 feet per 
second for the model would correspond to 100 feet per sec- 
ond for even a small airplane, such as the Weick Wl-A (ref- 
erence 2). It was thought that the discrepancy might be 
attributed to a gyroscopic torque (mentioned earlier) that 
was due to the interaction of the rotation of the wheel 
on its axle with lateral motion of the bottom of the tire. 
This lateral motion may be considered to be equivalent to 
a rotation, of part of the tire about an axis perpendicular 
to the spindle and wheel axle. This interaction gives 
rise to a spindle torque, which must be included in equa- 
tion (6). The numerical value of this torque will be 
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I „ 1 kp , where I ' is the moment of inertia about the wheel 

axle of the part of the tire that engages in the lateral 
motion (say, possibly l/4 of the total moment of inertia 
of the tire); is the angular rotation about the longi- 

* * 

tudinal axis = X/r); and q> is the angular- rotation 

about the axle (<p = V/r). Thus, the value of the gyro- 
static torque is I a XV/ r s ; but X = V ~ so, finally, 

the torque is I„ -^g The angular acceleration due to 

r ds 

J.-U- . Ip V 2 dX „, /-\ , ,,2 d 2 6 

this torque is ~~ — - — . Thus (6) 'becomes Y — - 

l v r s ds ds 

la 1 = o ^ _ - ^/ d-S Y Combining this equation with 

I w r 3 ds V Kj. / 

(9) 

If *!i + y s - ° 2 yS ** + JL^ ^ + -iS- l! ^ + 6 = 0 

Gx ds 3 MJi C x I"ir a i w Kj.' ds 3 I w r 2 C x ds 

Here the constants Gi, C2, and Ki are the same 
as before.' For the assembly on which the measurements 
plotted in figure 8 were made, I_ was 1.06 x 10** 4 slug- 
'feet , I & was 1.09 x 10~ 5 slug-feet 2 , and r was 

0.1S5 foot. Using these values, new theoretical diver- 
gences and frequencies were computed and these are com- 
pared with the experiment s in figure 10. It will be seen 
that, when the gyrostatic torques are included, the agree- 
ment of the theoretical divergences with experimental 
values is considerably improved at high velocities. 

The gyrostatic torque is proportional to dX/ds. 
Hence, at high velocities it is nearly proportional to 
d6/ds since -A and 6 are nearly in phase. Thus, this 
torque is approximately in the right phase to damp the os- 
cillation. 

Effect of Motion of Airplane 

When a shimmying wheel is attached to an airplane, 
the airplane will respond to a certain extent to the 
forces supplied by the shimmy. .This lateral motion of the 
nose of the airplane will in turn affect the shimmy and 
does, in fact, reduce the divergence at low velocities, 
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In order to calculate an approximation to this effect 
the airplane will he assumed to he a rigid body. By this 
assumption any deflection occurring in the structure or. at 
the rear tires is neglected. The gyrostatic torque due to 
lateral tire deflection is also omitted. 

The following symhols are used in the calculations: 

I z , moment of inertia about vertical axis through center 
of gravity of the airplane, 

m, mass of airplane. 

E;\ , lateral force exerted "by nose-wheel tire per foot de- 
flection. 

6, angle between central plane of wheel and average direc- 
tion of motion. 

\|/, angle between longitudinal axis of airplane and aver- 
age direction of motion. Measured positive in the 
same direction as 6. 

y, lateral distance between center of gravity of airplane 
and path of average motion. Positive in same direc- 
tion as X. 

Y , lateral force exerted by the rear wheels on the air- 
plane, 

l 2 and l lt horizontal distances from nose wheel and rear 
wheels, respectively, to center of gravity. 

1 = + X q . 

The side force developed by the shimmy is E^A. The 
lateral acceleration is 

m ds 

Similarly, the angular acceleration of the airplane is 
given by 

+ h Y = \[r = V S —~— (14) 
I z ds 2 
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There is, in addition, the geometrical relationship 

is. . ( e-*>b. (i 5 ) 

ds ds 1 I 

The lateral motion of the nose per unit forward motion, 
assuming that the rear wheels. do not skid, is 

ii + i — 

ds 

This motion will naturally affect the rate of deflection: 
of the nose wheel dA/ds; equation (9) "becomes 



dX c „ / 1 d6 \ . , d^ , 

ds 2 V K x ds/ ds 

Equation (6) is unaffected. Combining these equations 



+ T' k n ♦ h. V a 0.) + 1= 1 J £ 

i. K K X I / J ds 

+ (c x li 2 m 7 2 + I z V 2 Gx + I E A V 2 + I 8 l x ~) ~ " 

/ • I 

+ Ci flj m V 2 + -) 6=0 (17) 

\ Ki / 

This differential equation has two real roots, which 
are of little interest to the problem at hand, and two 
complex roots, which give the divergence and frequency of 
the shimmy under consideration. 

An experimental check on the theory was obtained "by 
running on a flat surface (see fig. 11) two models similar 
to the one sho\?n in figure 12. The model \vas launched "by 
a compressed-air catapult, which gave it a very smooth 
start. The nose wheel was locked during the acceleration 
so that the shimmy would not start to "build up before con- 
stant velocity was reached. The wheel was freed "by a lever 
tripped "by a wire. It then came into the field of view of 
an elevated camera and the flashlight "bulb made a track of 
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light on the camera film. The light entering the camera 
was interrupted 60 times a second "by a shutter driven by a 
synchronous motor. A typical trace is shown in figure 13. 
The curve was read only in the region "between the two ar- 
rows. Several such traces were made at various velocities 
and the divergence and frequency of the shimmy were meas- 
ured and plotted against velocity in figure 14. Two mod- . 
els with different mass characteristics and different caster 
angles and tires were used to get two checks on the theory. 
The theoretical divergences and frequencies in these models 
are also shown in figure 14. 

It will be noticed that both the theoretical and ex- 
perimental divergences are somewhat smaller at low veloc- 
ities. than in the case of the belt-machine shimmy. In 
general, however, the differences between the two cases arc 
small. The theory therefore seems adequate to account for 
the shimmy of freely rolling models when corrections are 
made for the response of the model to the shimmy. 

Pull— scale experiments similar to the ones with the 
model would be valuable in further checking the theory. 
One very rough check was made on the Wl-A. airplane (refer- 
ence 2). The friction in the spindle of this airplane was 
considerable but when the friction was reduced the wheel 
shimmied. Motion pictures were taken of the shimmying 
wheel and some of the enlarged frames are shown as figure 
15. The shimmy diverged and finally reached an equilib- 
rium amplitude. It will be seen that deflection of the 
bottom of the tire plays a prominent part in the shimmy. 
The frequency of the shimmy of this airplane was measured 
at a velocity of 47 feet per second and found to be 9 cy- 
cles per second. The frequency was then computed for this 
velocity (on the basis of constants found by the method to 
follow) and 10.5 cycles per second was obtained. These re- 
sults are taken to be indications that the theory will be, 
at least, approximately valid full scale. 

In the application of the theory, it is necessary to 
know the constants of the tire to be used. In the absence 
of such measurements for full-scale tires, it will be nec- 
essary to use the theoretical values of the parameters 

Ki = 2/r 2 and 0 2 = l/r. For I w C x (the spindle torque 
per unit lateral tire deflection for a completely asymmet- 
rical strain) the formula n r Ex will be used, where E\ 
is the lateral force per unit deflection, r is the radi- 
us of the tire, and n is an empirical nondimensional 



15 



factor to "be obtained from experiments. The factor n 
will, of course, increase with increasing caster length; 



ri for the model tires 
1.00 at caster angle 20 



was 
o 



0.69 at caster angle 5 and 



The 



se expressions for 0 if C s , and K x will Tbe used 



for numerical calculations in the rest of the report. 

METHODS OF AVOIDING- SHIMMY 
Spindle Damping 



In the past the only method that has been used to 
avoid shimmy in single swiveling wheels has consisted in 
the application of friction or hydraulic damping to the 
spindle. It will therefore be desirable to set up a cri- 
terion for the amount of spindle damping necessary to 
avoid shimmy. 

In this calculation and also in later calculations of 
other methods of avoiding shimmy the damping effect of gy~ 
rostatic torque due to lateral tire deflections will be 
neglected as well as that due to airplane motions; these 
omissions r.esult in small errors 'and usually give a con- 
servative results 

First, the effect of idealized viscous damping applied 
to the spindle will be calculated, that is, the effect of 
a torque pr.oportional to the angular velocity of the spin- 

die, say lit = KV^Ts" The angular acceleration resulting 
from this torque 7,7 --/ I w must be added to equation (6). 



-a d G 

v __ g 

ds 



_ r A l d6\ KV 



KV ae 
ds 



(18) 



Combining this equation with equation (9), which is unaf- 
fected by the addition of swindle friction, 



V a d 3 6 

Ci ds 3 VJx I 



) 



d 2 e 

ds 



K 



I 6 . + e = o 



I w ds 



(19) 



In the case of a self-excited oscillation represented 
by a third-order differential equation, Routh's discrimi- 
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nant provides a simple method of deciding whether the os- 
cillation will or will not diverge. If all the coeffi- 
cients arc positive (as they are here), then the oscilla- 
tion will diverge if and only if the product of the coeffi* 
cients of the third- and zero-order terms is greater than 
that of the first- and second-order terms. Hence for the 
shimmy to have zero divergence 



V s n s. 
o7 ~ \c7 il 



V /V K 1 Co 2 \ C 2 Z 



C x VJj I w K x G x / c x I w 
or, rearranging and solving for K/l w , 



iw Wiq 2 ; v Vv 5 ^! e x 3 / . c a 

for the condition of zero divergence. 

Only the positive value of K is of interest "because 
the application of negative damping to the spindle would 
.make the wheel unstable. From equation (20) it is found 
that there is a maximum positive value of K at the veloc- 
ity Jo 1 / C 2 2-l , which for the model wheels at 5° caster 
is 13.6 feet per second. Thus, it is seen (see fig. 9) . 
that the amount of spindle damping necessary to avoid 
shimmy is a maximum near the velocity at which maximum di- 
vergence occurs with an undamped wheel. This velocity 
will he denoted hy V Q . The damping required at this ve- 
locity, i.e., the maximum damping required, is 



Z max = W-T- L/ 1 + ^ 



CiCs / /., , Ki ,\ lh.-ft. 



- l) ^illTlll (21) 

' -rortinmc! n PT 55 ft C. . 



v v C 2 ' radians per sec. 



It will he seen from equation (21) that it is desira- 
hie to make the caster length small for in this way Ci , 
and hence K ma;x . , may ho made small. A rough approxima- 
tion to K na:x may he ohtained from equation (21) hy the 

use of the theoretical formulas for the tire parameters 
given earlier. Using those formulas and the value of n 
ohtained for 5° caster angle, 

W = 0.43 r 7 I w E A 



radians r>er sec. 
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Using these formulas, also 



Y 0 =.0.59 r s / -- ft. per sec. 

for the velocity at which this maximum damping is required. 
For the f 1-A airplane, using r = 0.67 ft., I w = 0.11 

slug-ft. a (estimated), and = 2,400 lb. /ft. (from 

rough measurements ) 

E nax = 4.6 lb. -ft. per radian per sec, and 
T 0 = 39 ft. per sec. 

For the Hammond Model Y airplane, which has had consider- 
able difficulty with shinny* r = 0.83 ft., I w (esti- 
mated) = 0.33 slug-f t. 2 , S x .(estimated) = 4,130 lb. /ft., 
which gave 

E„„^ = 13. Z lb. -ft. -ncr radian per sec, and 

T Q = 46 ft. per sec 

The amount of damping required increases with the 
fourth power of the size of the airplane. The application 
to an airplane of viscous spindle damping, such as might 
be applied by a hydraulic damper, night be expensive and 
rather heavy. It would bo necessary, moreover, to insure 
that adequate damping would bo obtained at all- amplitudes, 
that is, it would bo necessary to allow for the nonlinear- 
ity of the damper. 

Host practical shimmy dampers have used solid rather 
than viscous friction. Solid friction, however, is dif- 
ficult to deal with analytically in a straightforward man- 
ner. Therefore, resort will be had to the following sub- 
terfuge. It will be assumed that, if solid friction takes 
an amount of energy out of an oscillation equal to the en- 
ergy taken out by a given amount of viscous damping, it is 
equivalent to viscous damping for purposes of damping the 
oscillation. The energy taken out by viscous damping in 
l/4 cycle is 
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J nax 



J nax 



kg ae = 



kvH ae 
as 



whore G = G nax sin U) s. The energy taken out by a solia 
friction torque T is T 6 
the sane amount of energy 



max 



Therefore, to take out 



T = I KVU> G nax 



(22) 



(Telocity 
t an c e . ) 



enters "because the indepenaent variable is dis- 



Squation (22) gives an interesting explanation of the 
fact that shinny (and other self-excited vibrations) will 
usually not start in the presence of solid friction until 
an initial displacement is supplied. 1 Fron the theory of 
small oscillations the viscous friction required to damp a 
self— excited vibration is independent of the amplitude. 
(See equation (21).) It will he seen from (22) that the 
solid friction required to damp a self-excited vibration 
will he directly proportional to the initial amplitude of 
the oscillation. If solid friction is used to damp shimmy 
it will he desirable to insure that the wheel is centered 
when it first makes contact with the ground. This precau- 
tion has been found useful on airplanes. 



The value of or for the case of equilibrium viscous 
damping (zero divergence) is equal to the ratio of the co- 
efficient of the dG/ds term to the d 3 6/ds 3 term in 
equation (19). This fact may be readily verified by sub- 
stituting a simple harmonic solution for 8 and equating 
the sun of the inaginary terms to zero. Thus 



(JU 




(23) 



It is difficult to obtain a good estimate of the max- 
imum angle that should be used in computing the amount of 
solid friction necessary to damp shimmy on an actual air- 
plejie because the value depends on how much energy will be 
supplied to the shimmy after the tire starts to skid. If 
it be supposed that the energy supplied to the shimmy at 
angles greater than that at which skidding starts is en- 
tirely taken out by the solid friction, then the value of 
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6 for which skidding begins can "be put into (22) and a 
value of the friction will "be obtained such that shimmy 
with all initial angles will "be damped out. This supposi- 
tion is probably true for small caster lengths. The tire 
will start to skid when the lateral force is large enough 
to overcome the friction,, that is, when E^X = y,W, where 

|X is the coefficient of friction and W is the load on 
the swiveling wheel. Now, if the wheel moves in a direc- 
tion at an angle 9 with the central plane of the wheel 
the equilibrium deflection of the tire will be given by 

C 3 A = - 6 

from equation (8). In this way a rough estimate can be 
obtained of the angle 8 0 at which skidding starts. Thus 



G 2 u- W 

Substituting in (22) from (24) for 6 
(23) for uj, V n or 



max 



(24) 



the value from 



K 



max 



from (21) for 



yCi/CsKi for T , and the value of 
K, then 



T 

x nax 



-JLl_V_^ilL ( A + - if" (25) 



is obtained as the amount of solid spindle friction re- 
quired to damp shimmy at the velocity requiring most vis- 
■cous spindle damping (which velocity will probably re- 
quire the most solid spindle friction). The validity of 
formula (25) is restricted to small caster lengths; this 
restriction should not hamper its application because the 
caster length is preferably small. Using the theoretical 
formulas for the tire parameters, there is obtained the 
simple formula (again assuming small castor length), 

T nax = 0.1? n r I lb. -ft. (26) 

It should be remembered that the theoretical formulas 
for the tire parameters used in the derivation of (26) are 
approximately valid only for fires similar to the model 
tires used in the experiments. Thus (26) and all other 
formulas based on the theoretical tire parameters will be 
approximately valid only for air wheels or low-pressure 
tires. lor high-pressure tires it seems likely that the 
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effective r would be smaller than the geometrical radius 
of the tire. It should "be remembered that no conservative 
assumptions were made in the derivation of (26). 

If now |i = 0.55, a value of 19 foot-pounds spindle 
friction is required for the Wl-A airplane with W = 310 
pounds (including propeller thrust) and 44 foot-pounds 
for the Hammond Y airplane with W = 570 pounds. The 
amount of spindle friction actually present in the spindle 
and steering mechanism of the Wl-A was measured and found 
to be of the same order of magnitude as the computed value. 
This result explains why no shimmy was found with that air- 
plane. As mentioned previously, when this solid friction 
was reduced, shimmy occurred in the castering front wheel 
of this airplane. 

It will be seen that the amount of solid spindle fric- 
tion necessary to damp shimmy in nose wheels is considera-' 
ble even at small caster lengths. As is the case for the 
hydraulic damper, it increases with the fourth power of the 
size of the . airplane. 

Lateral Jreedom 

It has been seen that the divergence of the shimmy 
depends on the interaction of the lateral tire deflection 
with the angular rotation of the. spindle. If the coupling 
between these two degrees of freedom could be sufficiently 
reduced, the shimmy would not diverge. 

One method of reducing the coupling is to introduce 
a new degree of freedom into the oscillation. Thus, if 
the wheel is free to move in the direction of the axle, 
the tire deflections will depend on the inertia reactions 
of the wheel (and of any other members moving laterally)* 
The coupling present between A and 8 will depend on 
the mass of the wheel. The maximum mass that the wheel 
could have and still redtice the coupling sufficiently to 
avoid shimmy will now be calculated. 

If the lateral acc el erati on • of the wheel is x and 
the mass of the wheel and tire (and other members moving 
laterally) m', then 

E^A = m ' x = m 1 V 

,or 

/"» E^A a- . . 
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This lateral motion of the wheel must he added to 
equation (9 ) , 

dX , . / d6/ds\ dx 

• Ti =- e " ° 2 C x - ZTJ ' dl 

o r 

» = - 6 - C B ft - ^-^) - <H (28) 

ds 2 V E x / J m» V 

The effect on equation (6) will he small since only- 
small lateral motions will he involved. Comhining equa- 
tion (6) with equation (28), 

ll (V 2 ° 3 - + JL} ^ + + (i + J*) e = o (29) 

Oj. ds V Cj K x ' ds 2 Ciffl' ds v I^m'V ' 

It is seen that a damping term has heen added to the 
final equation of motion. (Cf. equation (10).) If this 
term is sufficiently large, it will counteract the diver- 
gence of the oscillation. Applying Eouth 1 s discriminant 
to ohtain the condition of zero divergence and simplify- 
ing, 

m 0 ' ~ ~ q-- slugs 

It is noted that m 0 ' is independent of velocity. If m 1 
is smaller than m 0 • the shimmy will converge; if great- 
er, it will diverge. Putting in the theoretical formulas 
for the constants 

m 0 1 = 1.45 -f slugs = 47 -~ pounds 

Jot the 171 -A airplane this value would he 11 pounds 
and for the Hammond Y , 22 pounds. Fortunately, these val- 
ues are of the same order of magnitude as the weights of 
the wheels used on these airplanes. 

The necessary amount of lateral freedom will now he 
estimated. This distance is determined hy the largest 
initial deflections that might he encountered. If the 
castor length is small, very little energy will he added 
to the shimmy after skidding starts, when IxX = \s. W. It 
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is necessary to find the amplitude of the oscillation in 
x for an oscillation in A with the amplitude 
Since m« x = E^X at maximum amplitude i 1 x = M.W. 

Suppose the value of m» used is that for equilib- 
rium damping (n 0 ' = CgE^/C^, then the maximum possible 
x is 



but 



h en c e 



X nax lA c 3 



x = V 2 tu 3 x 
•^max v w x nax 



x 



max a a _ 

As before, for equilibrium damping, uu 3 equals the ratio 
of the coefficient of the d8/ds term to that of the 
d 3 0/ds 3 term in equation (29), for m ' = Gq E^/Cj, , 



O) 3 = 



2 

C 2 T 



and 

x 



max 



Therefore, it will be seen that the minimum amount of lat- 
eral freedom necessary on either side of the wheel is 
equal to the _ amount that the tire will deflect laterally 
before skidding when the weight of the wheel is the re- 
quired value. This approximation is rough and no conserv- 
ative assumptions were made other than the one implicit in 
the neglect of whatever spindle friction may be present. 
For the Wl-A airplane the calculated amount of lateral 
freedom required is 0.85 inch on either side of the wheel 
and for the Hammond Y airplane, 0.91 inch. 

In the application of lateral freedom to an airplane 
wheel, means for keeping the wheel centered in the lateral 
travel must be provided. It can bo shown that the appli- 
cation of spring-restoring forces sufficiently large to 
accomplish this purpose would tend to make the wheel shim- 
my. It was found, however, in a lateral-freedom arrange- 
ment that, when the wheel was tilted, it would tend to 
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ride on the lowest portion of the axle. From figtire 16 
the point A has a smaller rolling radius than the point 
B but the point A must nake as many revolutions as 
point B. Hence, if no skidding occurs, point A will 
tend to lag behind point B . Then the resultant couple 
will turn the wheel about the spindle axis in such a direc- 
tion that the forward notion will sxipply a side force caus- 
ing the wheel to slide "do\7nhill" on the axle. This effect 
is found to be nore than large enough to overcome the com- 
ponent of the weight of the airplane that would tend to 
nake the wheel slide "uphill." Therefore, if an axle that 
is convex downward is used, the wheel will always seek the 
lowest point, which in this case will be the center of the 
lateral travel. Both theoretically and experimentally, 
this has been found to be a point of stable equilibrium. 

Two methods of applying lateral freedom to a cas-ter- 
ing wheel were tested on the belt machine. The first 
method was simply to allow the wheel to slide along the 
axle. The second method was to allow the wheel' to rotate 
about two axes in a vertical plane perpendicular to the 
axle. Shimmying in models could be avoided by either 
method. The second method, however, did not involve mak- 
ing the fork any wider. 

Additional full-scale tests of the first system were 
made of the Wl-A airplane. An axle curved downward was 
used to keep the wheel centered. A radius of curvature 
equal to 6 times the tire radius was arbitrarily chosen. 
It v/as found that the wheel centered satisfactorily and 
that on a rough turf-covered flying field shimmy was pre- 
vented by the use of l/2 inch of latteral freedom on ei- 
ther side of the wheel. On concrete with a lateral free- 
dom of 3/4 inch the wheel could be made to shimmy by taxy- 
ing over a severe bump at a moderate speed and an angle 
that would introduce a large initial , angular deflection of 
the wheel. In order to prevent shimmy under those condi- 
tions, the rotation of the swindle was limited by stops to 
±13°. The pilot did not regard this restriction of spin- 
dle freedom as objectionable for the ground handling of 
the airplane. 

The numerical values givon for the trials. of the Wl-A 
airplane are not to be taken as representative of the op- 
timum arrangement. The presentation of more conclusive 
design data must await the completion of further tests. 
One point worth noting is that recent model tests have in- 
dicated the possibility of dispensing with the use of 
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stops to limit the angular spindle notion, such as were 
found to he necessary in the trials with the Wl-A air- 
plane. 



CONCLUDING REMARKS 



The phenomenon called "kinematic shinny" seens to "be 
the fundamental phenomenon in the oscillation of nodel 
wheels. A dynamical theory "based on this discovery gives 
an adequate description of the shimmy of small model 
wheels. It seems unlikely that shimmy on full-scale air- 
planes will he fundamentally different. 

On the "basis of this theory a method has been pro- 
vided of numerically estimating the spindle friction nec- 
essary to damp shinny. An objection to this method" of 
avoiding shimmy lies in its indicated interference with 
the steering of nose wheels. 

A new method of avoiding shimmy by the use of lat- 
eral freedom for the wheel is introduced but requires fur- 
ther investigation before complete design data can be 
■given. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya. , August 11, 1937. 
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Figure 1.- Kinematic shimmy. 
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Figure 2.- Record of tbe track of a wheel executing a kinematic shimmy. 
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Figure 13.- Record from free-model apparatus. 



G3 



Direction of motion 



I II us t rated 
directions 
positive 



H Axle 



V- - Direc Hon o f mo tion 
Caster angle ^ 



Wheel central 
plane 



Spindle axis 



Axle 



Periph- 

era I 
center 

line 




Lateral displace- 
ment of typical - 
point = dX 




„ Wheel cenfral plane 



Path of typical point 
1 1 on tire relative to 

I 



wheel=ds 



Deformed peripheral 
/ center line 
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Figure 3.- Definition of terms. 
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Figure 5.- Geometrical relationships involved in kinematic 
shimmy; illustrating that d© = 2hh a <* Kl* 



Figure 6.- Illustrating the contribution 
of asymmetrical strain to dX/d 
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Figure 9.- Divergence and frequency of shimmy on belt 
machine compared with theory. 
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Figure 10.- Divergence and frequency of shimmy on belt 

machine compared with theory including gyrostatic 
torque. 




Figure 16.- Tilted 'axle. 
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Figure 11.- Apparatus for conducting free-madel experiments. (Elevation vis*), 



— i — i — i r~i i 

o Experimental a ivergence 




14 16 4 6 

Velocity, ft. per sec. 



8 10 12 

(a) Model 1 (*) 2 

Figure 14.- Divergence and frequency of shimmy of nose-wheel models compared with theox-y. 
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